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ABSTRACT. A method using a Green's function is developed for computing transient temperatures in 
a semitransparent layer by using the two-flux method coupled with the transient energy equation. Each 
boundary of the layer is exposed to a hot or cold radiative environment, and is heated or cooled by 
convection. The layer refractive index is larger than one, and the effect of internal reflections is included 
with the boundaries assumed diffuse. The analysis accounts for internal emission, absorption, heat 
conduction, and isotropic scattering. Spectrally dependent radiative properties are included, and transient 
results are given to illustrate two-band spectral behavior with optically thin and thick bands. Transient 
results using the present Green's function method are verified for a gray layer by comparison with a finite 
difference solution of the exact radiative transfer equations; excellent agreement is obtained. The present 
method requires only moderate computing times and incorporates isotropic scattering without additional 
complexity. Typical temperature distributions are given to illustrate application of the method by 
examining the effect of strong radiative heating on one side of a layer with convective cooling on the 
other side, and the interaction of strong convective heating with radiative cooling from the layer interior. 


1. NOMENCLATURE 

a absorption coefficient of semitransparent medium, m' 1 

c specific heat of semitransparent medium, J/kg-K 
Cj.-Cfi coefficients in Green's function 
D thickness of semitransparent layer, m 
coefficients in homogeneous solution 
F L blackbody energy fraction at t(X,r) in large frequency band 
F IL^2L blackbody fractions of external radiation fluxes in large frequency band 
g Green's function for radiative heat source 

G the flux quantity 2(q* + q;), W/m 2 * ; G = G/oT 4 
h,,h 2 convective heat transfer coefficients at x = 0, D, W/m 2 K 
Hi,H 2 dimensionless parameters, h,/oTf and h^oTf 
k thermal conductivity of layer, W/m-K 

m the quantity [3x^(1 - G)] 1/2 

n refractive index of radiating medium 

N conduction-radiation parameter, k/4oT-D 

P the quantity (2/3K D )[(l+p‘)/(l-p i )] 

q, radiative flux in the x-direction, W/m 2 ; q, = q/oT- 
q^,<£ radiative fluxes in positive and negative x directions, W/m 2 

external radiation fluxes oT^ and oT^ incident at x = 0 and D, W/m 2 
dimensionless radiation fluxes, q rl /oTf = (T.j/Tj) 4 and q r2 /oTf = (T./TJ 4 



T absolute temperature, K; t = T/T; 

Jgi.Ttf gas temperatures for convection at x = 0 and D, K 
t g „tj 2 dimensionless gas temperatures, T gl /Tj, 

Tj initial uniform temperature (used as a reference temperature), K 
T.i.T^ temperatures of blackbody radiative surroundings at x = 0 and D, K 
x coordinate in direction across layer, m; X = x/D 
0 time, s 

k d optical thickness of layer, (a + o,)D 
X c cutoff wavelength in vacuum, pm 

v radiation frequency, s' 1 

p density of semitransparent medium, kg/m 3 

p°,p‘ external and internal reflectivities at a boundary 
a Stefan-Boltzmann constant, W/m 2 K 4 

o, scattering coefficient in layer, m' 1 

x dimensionless time, (4oT-/pcD)0 

Q scattering albedo, o,/(a + a.) 

Subscripts 

L,S spectral bands with large and small frequencies 


2. INTRODUCTION 

For semitransparent materials where thermal radiation affects internal temperature distributions, transient 
behavior has been studied much less than steady state. There is, however, a significant need for detailed 
transient solutions to examine heat transfer and thermal stress behavior of ceramic components during 
reentry, thermal protection coatings, tempering of glass windows, porous burners and insulation systems, 
liquid drop space radiators, and some high temperature components in advanced aircraft engines. For 
these applications, internal emission is large, and unknown temperature distributions must be determined 
that are time accurate. The references considered here are for solid materials that have an index of 
refraction larger than one so internal reflections occur within the semitransparent region. 

An important early analysis was on the heat treatment of glass! Specular interface reflections were 
included, and two spectral absorption bands were used along with an opaque region for wavelengths 
larger than 4 pm. An analysis of transient temperature distributions in spacecraft windows for the 
thermal environment during reentry 2 followed somewhat along the earlier principles! A zonal method 
was used with a quadratic temperature distribution assumed in each zone. Some experimental 
verification was obtained by exposing single-pane windows to a high temperature radiation source and 
measuring transient surface temperatures. For a similar application, a recent study is for a silica window 
in high speed flow 3 Recent work 4 also provides additional experimental verification of transient 
temperature predictions for a glass plate. The plate was initially heated to a high temperature, and was 
then exposed to the ambient room on one side and to free convection and a radiant heater on the other 
side. Internal radiation effects were found to be very significant, as predicted by a two-band spectral 
calculation. Additional calculations for glass plates 5 used a nodal analysis based on a zonal method. 
General external boundary conditions were used with both radiation and convection. 

The analyses for glass plates have radiation combined with heat conduction, and scattering is negligible. 
In contrast a study with scattering that does not involve heat conduction is the cooling of a layer of liquid 
drops for dissipation of waste heat in the cold vacuum of outer space® Isotropic scattering was included 
within the droplet cloud. The resulting integral equation for the radiative source function was solved by 
iteration at each time step, and a forward finite-difference extrapolation was used in time. 

Many of the recent transient radiation analyses use some type of finite-difference method for solving the 
transient energy equation. The instantaneous local radiative heat source term in the energy equation has 



been evaluated by direct numerical integration’ by using the two-flux equations* by various expansion 
methods* 10 and with the differential approximation! 1 A semi-explicit finite-difference procedure with a 
■“fixed grid was used for a plane layer with scattering included! 2 Some effects of anisotropic scattering on 
transients have been examined 13 Transient solutions using the exact radiative transfer equations 7 are used 
to verify the accuracy of the present Green’s function method. These results are for a layer with 
refractive index larger than one with external convection and radiation at each boundary 7 

Most of the transient studies have been one-dimensional, usually for single plane layers. For a multilayer 
composite 1 , 4 a hybrid numerical method was used; the spatial variation in temperature was solved with 
finite differences, and a fourth-order Runge-Kutta method was used to extrapolate forward in time. For 
a two-dimensional rectangular region 1 , 1 a modified differential approximation was used for the radiative 
portion; a Crank-Nicolson finite-difference procedure was used for the transient energy equation. A fully 
implicit finite-difference solution 15 used an alternating direction implicit method for transient solutions 
in a two-dimensional square region. Two-dimensional Gaussian integration was used to evaluate the 
local radiative source in terms of the instantaneous temperature distribution. 

The present analysis develops a Green's function solution for the radiative source term using the two-flux 
equations. Various multi-flux methods have been discussed as a simplification for computing the 
radiative flux term in the energy equation 1 .* For the general boundary conditions of external convection 
and radiation on a layer with diffuse interfaces and unspecified surface temperatures, it was shown that 
the two-flux method can be used to predict accurate temperature distributions and heat fluxes for steady 
state! 7 By use of a shooting method to evaluate the radiative flux equations! the two-flux method was 
used to obtain transient solutions in materials with large refractive indices typical of ceramics. The 
shooting method, however, was limited to layer optical thicknesses less than about 5 where sufficient 
coupling exists between values at the two boundaries. Using the Green's function method, transient 
results are obtained here for optical thicknesses up to 40 and are compared with the opaque limit. An 
advantage of the two-flux method is that isotropic scattering is included without any additional 
complication. Typical transient results with large scattering are given to illustrate that scattering effects 
can be readily examined. The analytical relations are written for two spectral bands, and the form of the 
two-band equations shows how they can be expanded to include more detailed spectral variations. 


3. ANALYSIS 

A plane layer as shown in Fig. 1 has thickness D and radiative absorption and scattering properties that 
depend on frequency. The layer material is a semitransparent absorbing dielectric such as a ceramic, that 
is heat conducting and isotropically scattering with a refractive index n > 1. The layer boundaries are 
assumed sufficiently rough that they can be assumed diffuse!* For the results given here and for the 
nondimensionalization used, the layer is initially at uniform temperature T ; , but the theory applies for an 
arbitrary initial temperature distribution where T ; would then be a characteristic reference temperature. 
To begin the transient, the layer is subjected to surroundings so each boundary can receive radiative 
energy and is in contact with a gas that provides convective heating or cooling. Transient temperature 
distributions are obtained in the layer until steady-state is reached corresponding to the external radiation 
and convection conditions. 


3.1. Energy Equation and Two-Flux Relations for a Gray Layer 

The analysis is developed first for a gray layer and is then extended for spectral calculations. The energy 
equation in dimensionless form for the transient temperatures in the layer is 7 , 




Figure 1. Geometry and nomenclature for transient radiation and conduction in a semitransparent 

layer with spectral properties and isotropic scattering. 
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where -(l/4)aq r (X,t)/aX is the local radiative heat source arising from absorption, emission and 
scattering. Properties are assumed independent of temperature. The aq r (X,t)/aX is obtained in terms 
of t(X,t) from the two-flux relations which have been found to provide accurate results in previous 
steady-state 17 and transient® studies. The two-flux equations using the Milne-Eddington approximation 
are 1 , 6 


aq f (X,x) 

ax 


- k d ( 1- Q)[4nV(X,x) - G(X.t)] 


(2a) 


with <j(X,t) related to the radiative flux q^X,!) by, 


dG(X,x ) 

ax 


3 K D q r (X,t) 


(2b) 


The q, and Clare related to the positive and negative radiative fluxes shown in Fig. 1 by q, = - q£ and 

g=2(cj;+cd. 


3.2. Boundary and Initial Conditions 


For a semitransparent material there is radiant absorption only within the material interior so the 
conduction and convection conditions at the layer boundaries are 


at 

ax 


x*o 


- S [t *' ■ ,(0 - t)] 


at 

ax 


X=1 



[to,*) - y 


(3a,b) 


The radiative boundary conditions must include the effects of external and internal reflections at the 
surfaces that are subjected to dimensionless external diffuse fluxes q^ and q r2 . By considering the incident 



and reflected fluxes at and across an interface, the following boundary relations between G and q, were 
developed 1 , 7 

G(0,t) = 4 — q - 2 q (0,x) (4a) 

1-P‘ 1-P‘ 

G(1,t) = 4 — + 2 ^ q (l,x) (4b) 

1-p 1 17 l-p‘ 

To begin the transient solution, the specified initial condition for the results given here is a uniform 
temperature T(x,0) = Tj so t(X,0) = 1, but the method developed is valid for an arbitrary t(X,0). An 
initial distribution is also needed for dq/X.oydX. This is found by solving Eqs. (2a) and (2b) with 
t(X,x) = t(X,0), and the result is conveniently obtained from the Green's function solution that is 
developed. After substituting dqXXjOyaX into the energy equation, the transient solution is advanced 
one time increment by an implicit finite-difference solution. The new temperature distribution, t(X, Ax), 
is inserted into the Green's function solution to obtain dqXXjAxydX and the energy equation is then 
solved implicitly for t(X,2Ax). The process is repeated until steady conditions are achieved. The sections 
that follow will develop the Green's function solution and the implicit finite-difference solution of the 
energy equation. 


3.3. Green's Function for 3q r (X,x)/3X for a Gray Layer 

The Green's function is developed for a gray layer and is then extended in the next section for spectral 
calculations. Equation (2b) is differentiated in X and is used to eliminate 3es,(X,x)/dX from Eq. (2a); this 
yields a second order equation for G(X,x), 

~~Z ~ m 2 G(X,x) = -4m 2 nV(X,x) (5) 

dX 2 


where m 2 = 3 k d 2 (1 - Q). The x functional notation is omitted for convenience in what follows. The 
homogeneous part of Eq. (5) is satisfied by cosh(mX) and sinh(mX) so the Green's function must have 
the form 1 , 9 


g(X,0 = 


f C,(£)cosh mX + C 2 (£)sinh mX 
|c 3 (Ocosh mX + C 4 (Osinh mX 


0sx<5 

5<Xsl 


( 6 ) 


To obtain the coefficients C,(0...C 4 (0, the boundary conditions (4a,b) are written in terms of G by using 
Eq. (2b), 


G(0) = 4 


1-P° 

1-P 1 


+ 


2 l + p* dG | 
3x n 1- pi dX 


x=o 


(7a) 


G(l) = 4 


2 1 + p* dG 


1-p 1 3 k d i - p* dX 


(7b) 


x=i 


The g(X,Q in Eq. (6) must satisfy the homogeneous forms of Eqs. (7a,b), 


8 ( 0 ) 


- 2 l+P* ^ 
3k d 1- pi dX 


= 0 


g(l) + 


x=o 


2 1+ p* dg 

3k d l- p i dX 


= 0 


(8a, b) 


x=i 



( 9 ) 


By using Eq. (6) in Eqs. (8a, b) to evaluate C,(£)- ..C 4 (5X the g(X,£) 
C 5 (0 [ sinh mX + Pm cosh mX ] 

C 6 ($)[sinh m(l-X) + Pm cosh m(l-X)] 
where P = (2/3K D )[(l+p')/(l-p‘)]. 


g(x ,0 = 


is obtained in the form 
0iX<5 
5<Xsi 


The C } (5) and C 6 (0 are obtained by applying the two conditions 19 that g(X,£) is continuous at X = 5, 
and that with increasing X the derivative dg (X,£)/dX passes through a discontinuity of -1 at X = $ . This 
yields the two relations, 

C 6 (5)[sinh m(l-£) + Pm cosh m(l-l-)] = C s (£)[sinh m£ + Pm cosh m?] (10a) 


C 6 (5)[ _m cosh m(l -£) - Pm 2 sinh m(l-$)] 

- C s (5)[tn cosh mi; + Pm 2 sinh m£] = - 1 


(10b) 


Equations (10a,b) are solved for C s ($) and C s (5). and after substituting into Eq. (9) this gives the final 
form of the Green's function as, 


sinh m(l - £) + Pm cosh m(l-£) 


[m[(l +P 2 m 2 ) sinh m + 2Pm cosh m] 

x [ sinh mX + P m cosh mX ] 


0sX<5 


g(x,5) = 


(ii) 


sinh m£ + Pm cosh m£ 


m[ (1 + P 2 m 2 ) sinh m + 2Pm cosh m ] 


l <X <; 1 


*[ sinh m(l-X) + Pm cosh m(l-X)] 


The g(X,0 in Eq. (1 1) will be used to account for the nonhomogeneous term in Eq. (5) when computing 
G(X,t) at each time step during the transient calculations. To obtain the complete solution for G(X,t) 
the solution for the homogeneous part of Eq. (5) is needed. This is of the form G h (X) = E, sinh(mX) + 
E 2 cosh(mX). The boundary conditions Eq. (7a, b) are applied to evaluate E, and Ej with the result 

G h (X) = 4 l - {[sinh m(l -X) 

1- p‘ (1 +P 2 m 2 )sinh m + 2Pm cosh m (12) 

+ Pm cosh m(l-X)]q rl + [sinh mX + Pm cosh mXjq^J 


Then adding G h (X) and the nonhomogeneous solution obtained with g(X,$), the general solution of 
Eq. (5) at each time step is 

G(X,t)= G„(X) - 4m V [' g(X,Ot 4 (?,T)<l5 (13) 

Jo 

The (j(X,t) is substituted into Eq. (2a) along with t 4 (X,-c) to evaluate dq/X.tyaX which is used in 
Eq. (1) to solve for the temperature distribution at the next time step. 



3.4. Relations for a Two-Band Spectral Calculation 


" ^ relations developed for a gray material can be further developed for transient spectral solutions. To 
show how this is done, relations for two spectral bands are given. They are readily extended to more 
bands. Letting L and S designate bands with large and small frequencies, the radiative term needed for 
Eq. (1) is the sum of the large and small frequency contributions, 


aq r (X,T) _ aq rL (X,t) t dqJX.x) 

dx ex ex 

From Eq. (2a), 


aq fL (X,t) 

ax 


= *dlO- 0 L )[4nV(X,T)F L (t)- G L (X,t>] 


05) 


where F L (t) is the blackbody fraction at t(X,t) in the large frequency band. A similar relation is written 
for dqyax for the small frequency band. The Cj l (X,t) for Eq. (15) is found from Eq (13) written for 
the L band, 

Q l (X,x) = G hL (X) . 4m t Vy 1 g t (X,?)t 4 (?.t)F L (t)d5 (16) 

and similarly for <j s (X,t). The homogeneous solution needed for Eq. (16) for the L band is obtained 
fromEq. (12) as 


1 

(1 +PLniL)sinh m L + 2P L m L cosh m L 
x j[sinh m L (l-X) + P L m L cosh m L (l-X)]q rl F 1L 
+[sinh m L X + P L m L cosh m^Jq^F^} 


G hL = 4 


1-P° 

l-p 1 


(17) 


The Green's function in Eq. (16) for the L band is obtained from Eq. (1 1) as, 

sinh m L (l -£) + P L m L cosh m L (l -£) 


g L (XO = 


= < 


J^m L [(l +P L m L )sinh m L + 2P L m L cosh m L ] 

x[ sinh m L X + P L m L cosh m L X ] 


sinh m L \ + P L m L cosh m L \ 


[m L [(l +Pl m^sinh m L + 2P L m L cosh m L ] 

x [ sinh m L (l -X) + P L m L cosh m L (l -X) ] 


0 ±X<1 


(18) 


£<Xs 1 


Similar relations are written for the S band. With these relations Eq. (15) can be evaluated for the L 
band, and similarly for the S band; then dq.ldX is obtained from Eq. (14) for use in Eq. (1). 



3.5. Numerical Solution 


" Starting with the initial t(X,0), the Green's function is used to obtain G(X,0), and the dq,(X,0)/dX is 
evaluated from Eq. (2a) or Eq. (14) as described. Using Eq. (1), t(X,T) is integrated forward in time 
using the following implicit finite-difference algorithm at the interior points of a uniform grid: 


- — - T -t(X-AX,T + AT) + 
(AX) 2 


1 + 


« 

2NAt 

(AX) 2 


t(X,T + At) 


— t(X+AX,T+AT) 
(AX) 4 * 2 


- UX.T)- 


At d $ T 


4 9X 


X.T 


(19) 


At the boundaries Eq. (19) is applied and the t(-AX,x+At) and t(l+AX,T+At) that appear are 
eliminated by using the boundary conditions Eqs. (3a,b) in the forms. 


t(AX,T + At)- t(-AX,T+Ax) 
2 AX 
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4N 


L [t. 


- t(0,T + At)] 


(20a) 


«(1 AX,t * At) - ,(1 - AX,t + At) . . H, 

2 AX 4N ' 2 

This yields the special relations at X = 0 and 1 as, 


g2 J 


, 2NAt HjAt 
1 + + 


(AX) 2 2 AX 


2N At 

t(0,T + At) - — — — — t(AX,T +At) 
(AX) 2 
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4 ax 


I0,t 
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2AX 8 


(20b) 


(21a) 


_^H^l t (l-AX,T+At) + 
(AX) 2 


, 2NA t H 2 At 
1 + + 

(AX) 2 2AX 


t(l,T + At) 


= t(l,T)- 


Alfil + H 2 At t 

4 ax . 2 AX 82 

M* 


(21b) 


After advancing t(X,T) each At the radiant flux gradient is advanced to t + At by using the Green's 
function solution in Eq. (2a) or (14). The t(X,T) is then advanced to the next At using Eqs. (19) and 
(21). After checking various grid sizes it was found that 81 evenly spaced points across the layer gave 
accurate solutions for k d < 5. For k d > 5, 161 points were used which gave accurate solutions for k d s 
40 A time increment of At = 0.0025 provided stable and time accurate solutions for all results calculated 
here. 


4. RESULTS AND DISCUSSION 


The objective of this paper is to develop a convenient and accurate method for obtaining transient 

temperature distributions in a layer with isotropic scattering and spectrally dependent absorption and 

scattering properties. A Green's function solution was developed for the two-flux equations and this was 

used to obtain the radiative energy source term in the transient energy equation including externally 



incident radiation. The energy equation was solved with an implicit finite-difference method that included 
jconvective boundary conditions. The transient results that are presented serve two purposes. First the 
transient temperature distributions are compared with previous transient calculations using the exact 
equations of radiative transfer. This will verify the accuracy of the present transient two-flux Green's 
function method. The second puipose, as shown by the remaining figures, is to demonstrate some of the 
types of transient results that can be readily obtained with the present analytical and computational 
procedure. 

The transient temperature distributions given here begin with a uniform initial temperature t(X,0) = 1, 
although the solution applies for any t(X,0). Figure 2 shows typical comparisons of transient two-flux 
results for t(X,x) with those from a fully implicit method 7 using the exact transfer equations to evaluate 
the radiative flux gradient in the transient energy equation. The layer is heated at the hot side (X = 0) 
by a radiative flux equal to that from blackbody surroundings at T,, = 1.5Tj. At the other side (X = 1) 
the surroundings are at a low temperature T t2 = 0.5T ; so there is cooling by radiation through that 
boundary. The results in Fig. 2 illustrate the thermal behavior when there is convective cooling only on 
the side away from where radiative heating is applied. This simulates a possible condition for the wall 
of a combustion chamber when there is primarily radiative heating from combustion gases and soot on 
one side, and that side is not being film cooled. 

The three parts of Fig. 2 are for optical thicknesses k d = 2, 10 and 40 with no scattering and a refractive 
index, n = 2. The temperature distributions from the present analysis are solid lines; the dashed lines in 
Figs. 2a,b were calculated from the exact transfer equations! When x = 2 the t(X,x) are within 1% of 
steady state. The convective cooling at X = 1 produces a temperature decrease near that boundary after 
the transient begins. The results from the present method agree very well with those from the exact 
transfer equations; excellent agreement throughout the transient is obtained for k d = 10. 

The Green's function method can be applied for larger k d values than were used for the exact solution! 
Computational difficulties often arise for large k d because of the rapid local attenuation it provides. 
Transient temperatures for k d = 40 are in Fig. 2c. They are compared with t(X,x) for an opaque layer 
where radiation is absorbed and emitted only at the boundaries, and internally there is only heat 
conduction. These limiting results were obtained by removing the internal radiation terms from the 
computer program and adding radiative absorption and emission to Eqs. (3a,b) for opaque boundaries. 
For x = 2 the t(X,x) in the opaque layer is within 1% of steady state; the distribution is slightly curved 
and it becomes a straight line as x is further increased. The results show that for k d = 40 the layer is not 
thick enough to act as an opaque layer; there are very significant internal radiation effects on the transient 
and steady-state temperatures. A k d of at least 200 is required to approach the opaque limit. 

In Fig. 3 comparisons are made with the exact transfer equations for a very different set of boundary 
conditions. The conditions are the same on both sides of the layer so the t(X,x) are symmetric and only 
one-half of the distributions are shown. At the onset of the transient the gas temperatures are raised to 
two times the initial layer temperature, and there is strong convective heating at the boundaries with 
convection parameters H, = H 2 = 5. The surrounding radiative temperatures are low and, as the layer 
temperature increases, this produces radiative cooling from the layer interior. This causes the interior 
temperatures to remain much lower than the surface temperature even as steady state is approached. 
These results are for k d = 1 and they verify that the present method provides accurate solutions for a 
region of moderate optical thickness with large convective heating and interior radiative cooling. As 
shown in previous work 7 where results are given for n = 1, there are very significant effects of internal 
reflections when n = 2 as used here. 

An example of a transient two-band spectral calculation is in Fig. 4. At x = 0 a radiative flux of q,! = 3 4 
is applied at X = 0. The environment at X = 1 and the convening gases at both boundaries remain at the 
initial temperature, T^ = T^ = T^ = T ( . Hence, the layer is heated by radiation on one side with radiative 
cooling on the other side, and convective cooling at both sides. The steady state temperature distribution 



Dimensionless temperature, t(X,r) = T(X,r i )/T I 




(a) Comparison with exact results for k d = 2. (b) Comparison with exact results for k d = 10. 



(c) Comparison of k d = 40 with opaque layer. 

Figure 2. Two-flux, exact numerical, and opaque results for transient temperatures in a layer initially 
at uniform temperature after exposure to radiation on one side and convective cooling on the other. 
Parameters: Q = 0, N = 0.1, n = 2, q,, = 1.5 4 , q^ = 0.5 4 , H, = 0, H 2 = 1, tg = 0.5. 

for an opaque layer is the long-dashed line, and compared with this line the final profiles reached in the 
transient analysis show there is a strong effect of internal radiation. The short-dash and dot-dash lines 
are for gray layers with k d = 1 and 40. The solid lines are for two bands, k dl = 1 and k ds = 40 with A c T sl 











Dimensionless temperature, t(X,r) = T(X,r)/ T, 
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Dimensionless coordinate, X = x/D 

Figure 3. Two-flux and exact numerical results 
for transient temperatures in a layer initially at 
uniform temperature after exposure to hot gas 
convection on both sides with cold surroundings 
that produce internal radiative cooling. 
Parameters: Q = 0, k d = 1 , N = 0. 1 , n = 2, q rl = 
3i2 = 0.2 4 , Hj = H 2 — 5, t g i = tg2 — 2. 



Figure 4. Two-flux transient temperatures in a 
layer initially at uniform temperature after 
exposure to radiative heating on one side and 
convective cooling at both sides with T gl = T# = 
Tj; results for two bands compared with gray 
results for each band optical thickness. 
Parameters: Q = 0, N = 0. 1, n = 2, q,, = 3 4 , q r2 = 
1 , Hj = H 2 — 1 0, t gI = t g2 = 1 . 


= 4000 pm K; this band cutoff places about one-half of the incident radiation at T„ in each band. The 
general behavior of the two-band results is to be between the gray distributions at each time except for 
some deviation from this behavior at steady state. This demonstrates the importance of including spectral 
property variations. With the present method this adds very little complexity, and computer time 
increased about 50%. 


The two-band results are continued in Fig. 5 where scattering is introduced with Q = 0.9 or 0.99 in both 
bands. By raising £2 and keeping the same Kj* internal absorption is decreased and the transient 
temperatures do not respond as rapidly to the imposed radiative flux as for Q = 0. This also yields 
decreased temperatures at steady state. Since scattering is included in the two-flux relations, it does not 
add to the computer time or difficulty in evaluating the solution developed here. This is in contrast with 
needing to solve the scattering source function integral equation when using the exact transfer relations. 

In Fig. 6 the boundary conditions are of the type discussed for Fig. 3 where comparisons were made with 
results from the exact transfer equations. For all parts of Fig. 6 the layer is subjected for x > 0 to a high 
temperature gas at T^ = T^ = 4T ; so there is large convective heating on both sides. The surroundings 
are kept at T tl = T, 2 = Tj so radiative cooling occurs from the layer interior as the temperatures rise. 
Since the imposed conditions are symmetric, the t(X,r) are shown for one-half of the layer. The purpose 
of these figures is to illustrate that for a much different type of transient, spectral calculations including 
scattering can be readily carried out with the Green's function solution. The results in Fig. 6a are base 





Figure 5. Two-flux transient temperatures in a layer initially at uniform temperature after exposure to 
radiative heating on one side and convective cooling at on both sides with T gl ~ T g2 T;; results for two 
bands with large scattering compared with results for Q = 0. Parameters: N = 1, n = 2, cki = 3\ = 1, 

H, = H 2 = 10, t gl = 1,2=1, A C T„ = 4000 pm K. 

cases for gray layers with k d = 0. 1 and 10. The uniform steady-state temperature is also shown (dot-dash 
line) for an opaque layer (k d - «). For k d = 10 the temperatures at x = 1 (which are within less than 1% 
of steady state) are somewhat below the opaque limit except near the boundary; there is a small amount 
of radiative cooling from within the layer that provides lower interior temperatures. Radiative cooling 
has a very strong effect for k d = 0. 1 producing interior temperatures considerably below the surface 
temperature and the opaque limit. 

Retaining the same parameters. Fig. 6b demonstrates a two-band spectral calculation with k dl = 0.1, k ds 
= 10, and A c T,j = 4000 or 2000 pm K. Decreasing A c T fl places more energy in the k ds band, and the 
transient results change their behavior somewhat like increasing k d for the gray results in Fig. 6a. In Fig. 
6c scattering is included for Q = 0.9 and the parameters in Fig. 6b with A c T fl = 2000 pm K. Since 
increasing Q lowers the amount of absorption in the layer, there is an increased radiative cooling effect 
as was obtained for k d = 0. 1 in Fig. 6a. 


5. CONCLUSIONS 

A method was developed using a Green's function solution in conjunction with the two-flux equations 
to obtain transient temperatures in a semitransparent layer including internal radiation absorption, 
emission, and scattering with heat conduction. Starting from an initial temperature distribution that can 
have any arbitrary shape, the layer is heated or cooled by external radiation and convection. Transient 
temperature distributions from the present method were compared with results where the radiative flux 
gradient in the transient energy equation was evaluated numerically from the exact radiative transfer 
equations and a fully implicit finite difference solution was used. Good agreement was obtained for 




Dimensionless temperature, t(X,r) = T(X,r)/Tt 
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Figure 6. Two-flux transient temperatures in a layer initially at uniform temperature after exposure to 
convection by a hot gas on both sides, and cold surroundings that produce internal radiative cooling 
Parameters: N = 0. 1, n = 2, q,, = = 1, H, = H 2 = 5, t g j ^ t^ = 4. " 



optical thicknesses from 1 to 10 and for a refractive index of 2, including internal reflections from diffuse 
^boundaries. Two-band spectral calculations were carried out to demonstrate how the Green's function 
solution can be applied in spectral bands. The method can be readily extended to include more bands. 
Illustrative results also demonstrate the effect of isotropic scattering which is included in the present 
method without additional complexity. The technique is easy to program, has good numerical stability, 
and requires only modest computing times. It can be used to predict transient spectral behavior in 
materials subjected to large transient heating by radiation or convection starting from a specified arbitrary 
temperature distribution. 
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